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Abstract 

In the paper we discuss possible approaches to the problem of the rigorous derivation of quantum 
kinetic equations from underlying many-particle dynamics. For the description of a many-particle 
evolution we construct solutions of the Cauchy problems of the BBGKY hierarchy and the dual 
BBGKY hierarchy in suitable Banach spaces. In the framework of the conventional approach to 
the description of kinetic evolution the mean-field asymptotics of the quantum BBGKY hierarchy 
solution is constructed. We develop also alternative approaches. One method is based on the 
construction of the solution asymptotics of the initial-value problem of the quantum dual BBGKY 
hierarchy. One more approach is based on the generalized quantum kinetic equation that is a 
consequence of the equivalence of the Cauchy problems of such evolution equation and the BBGKY 
hierarchy with initial data determined by the one-particle density operator. 
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I. INTRODUCTION 

We develop a formalism suggested by Bogolyubov for the description of the evo- 

lution of infinitely many particles. The evolution equations of quantum many-particle sys- 
tems arise in many problems of modern statistical mechanics In the theory of such 
equations during the last decade, many new results have been obtained, in particular con- 
cerning the fundamental problem of the rigorous derivation of quantum kinetic equations 

221 1 and, among them, the kinetic equations describing the Bose 

n fi n 




condensate 



A description of quantum many-particle systems is formulated in terms of two sets of 
objects: observables and states. The mean value functional defines a duality between ob- 
servables and states. As a consequence, there exist two approaches to the description of the 
evolution. Usually, the evolution of many-particle systems is described in the framework 
of the evolution of states by the BBGKY hierarchy for marginal density operators 

HQ. 

An equivalent approach to the description of the evolution of many-particle sys- 
tems is given by the dual BBGKY hierarchy j^a] , [s(3] , j|J in the framework of the evolution 
of marginal observables. 

The aim of this work is to consider links between the many-particle quantum dynamics 
and quantum kinetic equations. 

A conventional approach to the problem of the rigorous derivation of kinetic equations 
from underlying many-particle dynamics consists in the construction of a suitable scaling 
limit jsj], for instance, the Boltzmann-Grad limit or the mean-field limit 13] , [2^] of a solution 
of the initial-value problem of the BBGKY hierarchy. As a result, the solution limit is 
governed by the limit hierarchy preserving the chaos property, and the one-particle density 



operator satisfies the kinetic equation 



1( 



Q, 



14j . In the paper we formulate new 



approaches to the solving the mentioned problem which are based on the description of a 
many-particle evolution by the dual BBGKY hierarchy and by the generalized quantum 



3 



kinetic equation. 

We outline the structure of the paper and the main results. 

In the beginning in Section II we introduce some preliminary definitions and construct a 
solution of the Cauchy problem to the dual BBGKY hierarchy for marginal observables and 
the canonical BBGKY hierarchy for marginal density operators of quantum many-particle 
systems. We formulate also one more approach to the description of quantum many-particle 
dynamics which is based on an equivalence of the Cauchy problem of the BBGKY hierarchy 
with initial data determined by the one-particle density operator and the corresponding 
initial value-problem for a generalized quantum kinetic equation. 

In Section III, the results obtained in the previous section are used to analyze of the mean- 
field asymptotics of constructed solutions, in particular to derive the nonlinear Schrodinger 
equation and its generalizations. We formulate also new methods of the derivation of quan- 
tum kinetic equations from underlying many-particle dynamics. One method is based on 
the study of the scaling limits of a solution of the initial- value problem of the dual BBGKY 
hierarchy. Another method is based on a generalized quantum kinetic equation. 

Finally in Section IV, we conclude with some observations and perspectives for the future 
research. 



II. DYNAMICS OF QUANTUM MANY-PARTICLE SYSTEMS 

We study possible approaches to the description of the evolution of quantum many- 
particle systems, namely the Heisenberg and Schrodinger pictures of the evolution. We 
introduce hierarchies of evolution equations for marginal observables and states and con- 
struct a solution of the Cauchy problems of these hierarchies in suitable Banach spaces. We 
develop also one more approach based on the generalized quantum kinetic equation that is 
a consequence of the equivalence of the Cauchy problems of such evolution equation and the 
BBGKY hierarchy for certain class of initial data. 



A. The dual BBGKY hierarchy 



We will consider a quantum system of a non- fixed (i.e. arbitrary but finite 23|) number 
of identical (spinless) particles obeying the Maxwell-Boltzmann statistics in the space M. u . 
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We will use units where h = 2irh = 1 is the Planck constant and m — 1 is the mass of 
particles. The Hamiltonian of such a system H = @^ =a H n is a self-adjoint operator with 
the domain V(H) = = ©~ ^ n G T H | ^ € V(H n ) G H n , J2 n \\ H n^n\\ 2 < 00} C .F w , 
where JF W = 0^0?^®™ is the Fock space over the Hilbert space Tt (Ti? = C). Assume 
Ji = L 2 (W U ) (the coordinate representation), then an element G T-u = Q)^ =0 L 2 (W n ) 
is a sequence of functions ip = (ipo, ipi(qi), ■ ■ ■ , VVi(<Zij • • • , Qn), ■ ■ ■ ) such that \\ip\\ 2 = 
l^o| 2 + Y^=i I dqi ■ ■ ■ dq n \ifj n (qi, . . . , q n )\ 2 < +00. On the subspace of infinitely differen- 
tiable functions with compact supports i\) n G L\{W n ) C L 2 (W n ) the n-particle Hamiltonian 
H n acts according to the formula (H = 0) 

n n 

H n ip n = K( S)^n + € ^ (1) 

1=1 i<j = l 

where K(i)if) n = —^A q .ip n is the operator of kinetic energy, ^{i,j)tp n = — qj\)ip n is 

the operator of a two-body interaction potential satisfying Kato conditions, and e > is a 
scaling parameter. 

Let a sequence g = (<7o> <7i> • • • ? <7tu • • • ) be an infinite sequence of self-adjoint bounded 
operators g n defined on the Fock space Tu = 0^L o 7i® n over the Hilbert space H (H° = C). 
An operator g n , defined on the n-particle Hilbert space H n = H® n will be denoted by 
g n (l, . . . ,n). For a system of identical particles obeying the Maxwell-Boltzmann statistics, 
one has g n {l, . . . , n) = g n (ii, ■ ■ ■ ,i n ) f° r an Y permutation of indices . . . , i n } G {1, . . . , n}. 

Let the space 2(J 7 n) be the space of sequences g = (g ,g 1 ,..., g n , . . . ) of bounded 
operators g n defined on the Hilbert space Ti n and satisfying the symmetry property 
g n (l, . . . , n) = g n (ii, . . . ,i n ), if . . . ,i n } G {1, . . . , n}, with an operator norm. We will 
also consider a more general space £ 7 (JF W ) with a norm 

IMk(^) = ™f x — s ll^nlkWn), 

where < 7 < 1 and ||.||,c(ft„) is an operator norm. An observable of the many-particle 
quantum system is a sequence of self-adjoint operators from £ 7 (jF w ). 

On the space £ 7 (JF W ) we consider the initial- value problem of the dual BBGKY hierarchy. 

The evolution of marginal observables is described by the initial-value problem of the 
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following hierarchy of evolution equations: 

a s s 

Q-G s {t,Y) = (^AT (z)+ e J2 Mi*(i,j))G.(t,Y) + 

i=l i<j=l 

s 

+e E ^UJuh)G s ^(t,Y\{j 1 }), (2) 

G s (t) |t= = ^(0), s > 1. (3) 

In equations fl2} we use notation Y = (1, . . . , s). The operators Ao, J\f int are consequently 
defined on X>(Ao) C £ 7 (J-^) as follows: 

^oCy)» = -i[^^(«)], (4) 
M a t(i,j)g = -i[g,&(i,j)], (5) 

where [ • , • ] is a commutator of operators. We refer to the evolution equations (j2J) as the 
quantum dual BBGKY hierarchy, since the canonical BBGKY hierarchy [5] for marginal 
density operators F(t) is the dual hierarchy of evolution equations with respect to the 



following bilinear 
observables) 26| , 



orm - average values of observables (mean values or expectation values of 



3l|: 



oo 

(G(t), F(0)) = £ -Tri,..., s G s (t)F s (0). (6) 



si 

s=0 



If TC — L 2 {R V ), the evolution equations (|2J) in terms of kernels of the operators G s (t), 
s > 1, are given in the form of the equations 

d 

i-^G s (t,q 1 ,... : q s ;q 1 ,...,q s ) = 

^ s s 

= ( - 2 X^ _A * + A <^ + e S ~ $ ~ _ ?i, si, ei) + 

1=1 l=i<jr 

s 

+e E - gj) - - qj ))G s ^{t, q 1} . . . , V, . . . , g s ; q[, . . . , V, . . . , q' s ), 

j 

where (q u . . . , V, . . . , q s ) = (q u . . . , q^i, q j+1 , ...,&). 

To construct a solution of the abstract initial- value problem (J2J)- (J3|) we introduce some 
necessary facts. 
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If g G we introduce the group Q{t) = ®^ =0 G n (t) of operators 

Q n (t)g n =e UH "g n e- ltH ". (7) 

A solution of the initial-value problem of the Heisenberg equation for observables of quantum 
many-particle systems is defined by this group of operators. 

On the space £ 7 (JF W ), the one-parameter mapping M} 3 t \— > Q(t)g defines an isometric 
*-weak continuous group of operators, i.e. it is a Cg-group. The infinitesimal generator 
Xf = ©^L Xf n of the group of operators (JTj) is a closed operator for the *-weak topology 
and on its domain of the definition T>(N) C £ 7 (JF W ), which is everywhere dense for the 
*-weak topology, Kf is defined in the sense of the *-weak convergence of the space £ 7 (.Fft) 
as follows: 

w*-limUg(t)g-g)=i(Hg-gH)=Mg. (8) 

Here, H = ($^ = qH ti is Hamiltonian (Tj[|) of the many-particle system and the operator: 
Mg = —i(gH — Hg) is defined on the domain T>(H) C T-u- We remark that operator (jSJ) is 
the generator of the Heisenberg equation. 

We introduce the nth-order (n > 1) cumulant of the groups of operators d7j) [27J , [28j , [29j 

X n {t)=f* n (t,X)= J2 (-1) |PM (|P|-1)! II ( 9 ) 

where is the sum over all possible partitions P of the set X = (l,...,n) into |P| 
nonempty mutually disjoint subsets Xj C X. 

We indicate some properties of cumulants (jSJ) of groups of operators (JZJ) 32]. If n = 1, for 
gi G T>{M\) C £(7ii), the generator of the first-order cumulant in the sense of the *-weak 
convergence of the space £(Hi) is given by operator (jSJ), i.e. 

w*-lim (fai{t, 1) - I)gi - (Afg)^ = 0, 

where the operator Xf is defined by (151) or (jl)). In the case n = 2 we have, in the sense of 
the *-weak convergence of the space £(7^2), 

w*- lim (- 2t 2 (t, 1, 2)g 2 - 6(^(1, 2) )g 2 ) = 0. 

If n > 2, as a consequence that we consider a system of particles interacting by a two-body 
potential ([I]), it holds 

w*-limi2l n (t)5f n = 0. 



We introduce also some abridged notations: Y = (1, . . . , s), X = Y\{ji, . . . , j s -n}, the 
set (Y\X)i consists of one element of Y\X = (ji, . . . ,j s - n ), i-e. the set {ji, . . . ,j s - n } is a 
connected subset of the partition P (|P| = 1, |P| denotes the number of partitions). We will 
also denote the set (Y\X)i by the symbol {ji, . . . ,j s - n }i- 

On the space £ 7 (JF W ) for the abstract initial- value problem (J2])-(J2D the following statement 
is valid. 

A solution of the initial-value problem to the quantum dual BBGKY hierarchy (d))-([3|) is 
determined by the expansion (s > 1) 

s 1 s 

G s (t,Y) = Y,T 7 -_ y E %i+n{t,(Y\X) 1 ,X)G s ^ n (0,Y\X), (10) 

where the operator 2l 1+n (t, (Y\X) 1 , X^j is the (1 + n)th-order cumulant (TJ|) defined by the 
formula 

z 1+n {t t (y\x) lt x) = (-i) |PM (|P|-i)! n 

p : {(y\x) 1 ,x}=UiX !i x,cp 
// G(0) G T>{N') C £ 7 (JF W ) ; it is a classical solution, and, for arbitrary initial data 

G(0) G £ 7 (JF W ) ; it is a generalized solution. 

Thus, solutions of the first two equations of hierarchy ([2]) are given by the following 

expansions 

G 1 (t,l)=% 1 (t,l)G 1 (0,l), 

G 2 (t, 1, 2) = Sti (t, {1, 2}i)G 2 (0, 1, 2) + 2l 2 (t, 1, 2) (^(0, 1) + G?i(0, 2)) , 
where the first-order cumulant 2ti (t, {1, 2}i) = ^2^, 1, 2) is defined by group ([7j). 

B. The BBGKY hierarchy 

The sequence F — (I, Fx, ... , F n , . . .) defined on the Fock space JF W of self-adjoint positive 
density operators F n (I is an identity operator) describes the state of a quantum system 
of a non-fixed number of particles. The marginal density operators F n , n > 1, whose 
kernels are known as marginal or n-particle density matrices, defined on the n-particle 
Hilbert space 7i n = 7i® n , we denote by F n (l, . . . ,n). For a system of identical particles 
described by the Maxwell-Boltzmann statistics, one has F n (l, . . . ,n) = F n (ii, . . . , i n ) if 
Oi, • -.,4} G {l,...,n}. 
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We will consider states of a system that belong to the space S} a {T-H) = ©^Lo^" ^ C%i) 
of sequences / = (l,fx,...,f n ,...) of trace class operators /„ = f n (l, . . . , n) e £ 1 ('H n ) 
satisfying the above-mentioned symmetry condition, equipped with the trace norm 

oo 

ll/lk(^) = E« nTr wl/n(l,...,n)|, 

n=0 

where Tri r .. jn is the partial traces over 1, . . . , n particles, and a > 1 is a real number. By , 
we denote the everywhere dense set in S^ Jl {J-'u) of finite sequences of degenerate operators 
with infinitely different iable kernels and compact supports. 

On the space we consider the following initial- value problem of the quantum 

BBGKY hierarchy 



^-F s (t) = -(J2^ + e E ^(i,j))F s (t) + 

i=l i<j=l 

s 

+ Tr s +i( - s + l))F s+1 (t), (11) 

i=l 

F s {t) \ t=0 = F s (0), s > 1. (12) 

If / G £q(JF w ) c V(Af) C £^(JF W ), the operators Mq, Afmt ar e consequently defined by (j4j), 
([H]). We remark that hierarchy (fTTj) is the dual hierarchy of equations to hierarchy ([2]). 

In terms of the kernels F s {t, (ft, . . . , q s ; q' 1: . . . , q' s ) of s-particle density operators F s (t), 
i.e. marginal or s-particle density matrices, equations (llip take the canonical form of the 
quantum BBGKY hierarchy 



. d _ . / /. 

^F s {t,q 1 ,...,q 3 ;q 1 ,...,q s ) 



+ / rf ^+i(^(% - g, + i))-Fi+i(*, 9i, g s +i; gl,..., i,? a +0- 

i=l ^ 

To construct a solution of the initial- value problem (fTT|) - ([T2"|) . we introduce some prelim- 
inary facts. 

On the space £^(JF W ) we define the following group G(—t) = ®^=oGn(—t) of operators: 

Sn{-t)f n :=e- UH " f n e UH ". (13) 



On the space £^(JF W ), mapping (|T3|) : t — > Q{—t)f is an isometric strongly continuous group 
which preserves the positivity and the self-adjointness of operators. A solution of the initial- 
value problem of the von Neumann equation for a statistical operator is defined by this 
group. 

If / G ^aoi^Fn) C V(Af) in the sense of the norm convergence of the space 
there exists a limit, by which the infinitesimal generator —A/" = ©^L (— A/" n ) of the group of 
operators f|T3|) is determined as 

Jim ~ (G(-t)f -f)= -i(Hf - fH) := -A/7, (14) 

where H = ®°^ =0 H n is Hamiltonian ([1]) and the operator —i(Hf — fH) is defined on the 
domain T>(H) C Tu- We note that operator (THl) is the generator of the von Neumann 
evolution equation. 



Let X = (1, . . . , n). The nth-order cumulant [28(,[29|] of the groups of operators ( 1131) is 
defined as (n > 1) 

f*n(-t)=f& n (-t,X)= J2 (-1) |PM (|P|-1)! II G\ Xi \(-t), (15) 

P:X=U l X l X t CP 

where Y2p is the sum over all possible partitions P of the set {1, . . . , n} into |P| nonempty 
mutually disjoint subsets X$ C X. 

If n — 1, for fx e £j(7ii) C V(Mi) C £ 1 (7ii) in the sense of the norm convergence in 
£ 1 (?ii), the generator of the first-order cumulant is given by operator (tHj), i.e. 

SS U^W-*, ^ " J )* " ( - ^)lL(Wx) = °- 
In the case n = 2, for cumulant ffT5l) we have 



lim |||saa(-*)/ a - e( - Ai nt (l, 2))h\\ &m = 0. 
If n > 2, as a consequence that we consider a system of particles interacting by a two-body 



potential, it holds 



!™II^H)^Li(^) = - 



We introduce the following notations: Yp = (X\, . . . ,X\p\) is a set, whose elements are 



|P| mutually disjoint subsets Xi C Y = (1, . . . , s) of the partition P : Y = Ui=i-^i- Since 
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Yp = (X±, . . . ,X\p\) then Y\ is the set consisting of one element Y = (1, ...,s) of the 
partition P (|P| = 1). To underline that the set (1, . . . , s) is connected subset (the cluster 
of s elements) of a partition P (|P| = 1), we will also denote the set Y\ by the symbol 

On the space £} a (Tu) for the abstract initial- value problem (TTTT) - (HT^T) the following state- 



ment is valid 



29]. 



If F(0) G £*(JF W ) and a > e, then, for t G M. 1 , there exists a unique solution of initial- 
value problem < f77]) -< f7^) given by the expansion (s>l) 

F.{t,Y) = ^-Tr s+1 ,... )S+n 2l 1+n (-t,y 1 ,s + l,...,s + n)F s+n (0,X), (16) 

where 



n=0 



Sti +n (-*,yi,a + l,...,a + n)= £ _ i)i JJ £ |Xi| (-t) 

P:{Yi,X\y}=U i X i X lC P 

is the (1 + n)th-order cumulant < T73]) o/ £/ie groups of operators l[TS\) . Y^ P is i/ie sto over 
a// possible partitions P o/ i/je set {Yi, s + 1, . . . , s + n} into |P| nonempty mutually disjoint 
subsets XiC {Y U X\Y}. 

For initial data F(0) G £^ (JF W ) ; it is a strong solution, and, for arbitrary initial data of 
the space £^(JF W ) it is a weak solution. 

The condition a > e guarantees the convergence of series (fTBI) and implies that the mean 
value of a number of particles is finite. This fact follows if we renormalize sequence ffTB]) in 
such a way: F s (t) = (N) s F s (t). For arbitrary F(0) G ^(JF-^), the mean value (jEJ) of the 
number of particles 

(7V)(t) = Tr 1 F 1 (t,l) (17) 

in state ([TBI) is finite. In fact, 

IW(*)l<CaimO)lk(^)<00, 

where c a = e 2 (l — ^) _1 is a constant. To describe the evolution of an infinite-particle system, 
we have to construct a solution of the initial- value problem ([TTi) - (lT2l) in more general spaces 
than £^,(JF W ). This problem will be discussed in Conclusion. 

We remark that, for classical systems of particles, the first few terms of cumulant expan- 
sion ([TBI) for the BBGKY hierarchy were obtained in [z|],[sj]. The methods used by Green and 
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Cohen were based on the analogy with the Ursell-Mayer cluster expansions for equilibrium 
states. 



A solution of the initia 
(iteration) series 



MM 



value problem (TTT1) - (IT2|) is usually represented as the perturbation 



12J. On the space ^(J 7 ^), expansion ffTol) is equivalent to the 
iteration series. Indeed, if an interaction potential is a bounded operator, then if f s G £ 1 (7Y S ), 
for group (fl3|) . an analog of the Duhamel formula holds 

{g 8 (-t,i,...,8)-]lg 1 (-t,i))f a = 

1=1 

Is s 

= e / drl[g 1 (-t + T,l)(- J2 ^Uh3))Qs{-r)f s . (18) 
o *=i »<i=i 

Then, according to the unitary property of group f|T3|) on the space £q(JF w ), solution ex- 
pansion (Tl6l) reduces to the iteration series of BBGKY hierarchy (TTTT) 

i *n-l 

F 8 (t, 1, . . . , s) = J] e n / dh . . . I db n Tr s+lt .„ !S+n g s {-t + h) x 

n=0 ^ ^ 

s 

^2 (-Nmtih, S + l))G s +l(-tl +t 2 ) ... Gs+n-l{-tn-l + t n ) X 



X 

h=l 

s+n—1 

X 

i„ = l 



(- Mnt(i n ,s + n))G s+n (-t n )F s+n (0,l, . . . ,s + n). (19) 



If F(0) G £q(JF w ), this series exists and converges for a finite time interval \&\, 11]. 

As was mentioned above, functional ([6]) of mean values defines a duality between marginal 
observables and mare; inal states. If G(t) G £ 7 (.F W ) and F(0) G £«(^w), then functional © 
exists, provided that a = 7 _1 > e, and the following estimate holds: 

\(G(0),F(t))\ = \(G(t),F(0))\ <e 2 (l- 7 c)- 1 ||C?(0)|| Jly( ^ ) ||F(0)|| i j i _ i( ^ ) . 
C. The generalized quantum kinetic equation 

We consider one more approach to the description of the evolution of states of quantum 
many-particle systems. Let the initial data are completely characterized by the one-particle 
density operator iq(0), for example, the initial data satisfy the chaos property (Maxwell- 
Boltzmann statistics) 

fM(o) = (/,iq(o,i),. ..,n* = /i(o, *),-.. )• 

12 



In that case, the initial-value problem of BBGKY hierarchy (fTTl) - (fT2"]) is not a completely 
well-defined Cauchy problem, because the generic initial data are not independent for every 
density operator F a (t), s > 1, of hierarchy of equations (TTTT) . Thus, it naturally yields the 
opportunity of reformulating such initial-value problem as a new Cauchy problem for the 
one-particle density operator, i.e. Fi(t), with independent initial data i*\(0) and explicitly 
defined functionals F s (t, 1, . . . , s \ Fi(t)) , s > 2, of the solution Fi(t) of this Cauchy problem 
instead other s-particle density operators F s (t), s > 2 j^.j^. 

Consequently, for an initial state satisfying the chaos property, i.e. F( c )(0), the state of 
a many-particle system described by the sequence F{t) = (I, Fi(t, 1), . . . , F s (t, 1, . . . , s), . . .) 
of s-particle density operators (TTB]) can be described by the sequence 

F(t | F 1 (t)) = (I, F x (t, l),F 2 (t, 1, 2 | F 1 (t)) ,...,F s (t,l,...,s\ F x (t)), . . . ) 

of the functionals stated above. 

At first, we define the sequence F(t \ Fi(t)) of functionals. The functionals F s (t, 1, . . . , s \ 
Fi(t)), s > 2, are represented by the expansions over products of the one-particle density 
operator Fi(t) (for particles obeying Maxwell-Boltzmann statistics) 

oo ^ s+n 

F a (t,l,...,s | Fi(t)) =^-Tr s+1 ,... iS+n Q3 1+n (t) J]Fx(t,i), (20) 

n=0 ' i=l 

where the evolution operators 9Ji + „(t) = 23i +n (t, {1, . . . , s}i, s + 1, . . . , s + n), n > 0, are 
defined from the condition that expansion (1201) of the functional F s (t \ F 1 (t)) must be equal 
term by term to expansion (fl6|) of the s-particle density operator F s (t). 
The low-order evolution operators 23i +n (t), n > 0, have the form 

aJi(t,yi) = a 1 (t,ro J (21) 

s 

QJ 2 (t, Y 1 ,s + l)=% 2 (t,Y 1 ,s + l)- Mt, Y 1 ) Mt,j, s + 1), (22) 
where 2t n (t) is the nth-order cumulant (semiinvariants) of scattering operators 

n 

g n {t,l,...,n) :=g n (-t,l,...,n)Y[g 1 (t,i), (23) 

i=i 

Qi{t) = I is the identity operator. 



13 



In terms of scattering operators (1231) evolution operators (1211) . (122"]) get the form 
Wi(t,Y 1 )=G a (t,Y), 

s 

oj 2 (t, y l5 s + 1) = g s+1 (t, Y,s + i)- g s (t, y) £ j, s + 1) + (* - 1) r). 

i=i 

For iq(0) G ^(W) the sequence F(t | Fi(t)) of functional (J2DJ) exists and series fl20|) 
converges under the condition that ||Fi(0)|| < e _1 , i.e. if the mean value of particles is finite 



25]. 



We remark that expansions (I2(jp are an nonequilibrium analog of expansions in powers 
of the density of the equilibrium marginal density operators (^{Jls]. 

We now formulate the evolution equation for the one-particle density operator Fi(t), i.e. 
for the first element of the sequence F(t \ Fi(t)). If ||Fi(0)|| < e _1 , it represents by series 
ffTB"]) convergent in the norm of the space S}(7i) 

oo ^ n+1 

Fi ^' v = J2 d Tr2 -" i+n i; •••'"+ v n «)» ( 24 ) 

n=0 ' i=l 

where 2ti +n (— £) is the (1 + n)th-order cumulant (|T5|) of groups of operators (fT3l . Let 
Fi(0) G £,l(7i). Then, by differentiating series (|24p with respect to time variable in the sense 
of the norm convergence of the space SI 1 (Hi), according to properties of cumulants (IT5|) . we 
find that the one-particle density operator Fi(t) is governed by the initial- value problem of 
the following nonlinear evolution equation (the generalized quantum kinetic equation) 

-Fi(t : 1) = -Afi(l)Fi(t, l) + Y, ^ Tr 2,3,..,n+2 ( - Mnt(l, 2))3J 1+n (t) J] F 1 (t, z), (25) 

n=0 ' i=l 

F 1 (t,l)| t=o = F 1 (0,l). (26) 

In the kinetic equation (125]) . the evolution operators QJi +n (t) = QJi +n (£, {1, 2}!, 3, . . . , 2 + 
n), n > 0, are defined as above. 



For initial- value problem (|2"3"l) - (j2"6"l) the following statement holds 25]. 
If Fi(0) G 2,q(H) is a non-negative density operator, then, provided ||Fi(0)|| < e _1 , there 
exists a unique strong global in time solution of the initial-value problem ( d5j) -( PS| ) which 
is a non-negative density operator represented by series [2$) convergent in the norm of the 
space £}(1-L) and a weak one for arbitrary initial data .Fi(O) G £ 1 (7^). 
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As a result, the following principle of equivalence of the initial- value problems (TTTT ) - (TT2|) 
and (|25]-(P5|) is true. 

If the initial data are completely defined by the trace class operators Fi(0), then the Cauchy 
problem / f77]) -/ f7^j) is equivalent to the initial-value problem I[2d*\) - I[2o}) for the generalized 
kinetic equation and functionals F s (t, l,...,s | F 1 (t)), s > 2, defined by expansions / TUZTj) 
under the condition that ||Fi(0)|| < e _1 . 

We note that this statement is valid also for more general initial data than F^(Q), 
namely the initial data determined by the one-particle density operator F\ (0) and operators 
describing initial correlations. In this case the initial correlations are a part of the coefficients 
of equation (I2"5"j) and functionals (|2"Uj) . 

Thus, if the initial data are completely defined by the one-particle density operator, then 
all possible states of infinite-particle systems at an arbitrary moment of time can be described 
within the framework of the one-particle density operator without any approximations. 

We remark that functionals ( |20l) are formally concerned with the corresponding function- 
als of the Bogolyubov method of the derivation of kinetic equations [1] . Indeed, functionals 
( |20l) and corresponding Bogolyubov functionals coincide if the principle of weakening of cor- 
relations for functionals (1201) holds. The proof of this assertion is similar to the proof jlj] of 
an equivalence of the BBGKY hierarchy solution (jT6l) and iteration series (|T9l) . 

III. DERIVATION OF NONLINEAR SCHRODINGER EQUATION 

We consider the problem of the rigorous derivation of quantum kinetic equations from 
underlying many-particle dynamics by the example of the mean-field asymptotics of above- 
constructed solutions of quantum evolution equations. In subsections B and C we formulate 
new approaches to the derivation of a nonlinear Schrodinger equation. 

A. The mean-field limit of the BBGKY hierarchy solution 

We present the main steps of the construction of the mean-field asymptotics of solution 
( I16p of the initial- value problem (I11I) - (I12I) . For that, we introduce some preliminary facts 
on the asymptotic perturbation of cumulants. 

If f s G £ 1 (7i s ), then for an arbitrary finite time interval, there exists the following limit 
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of strongly continuous group (flcfl) : 

s 

lim|| (&(-*) -tlGi{-t,j))f.\\#M = 0. (27) 
i=i 

According to an analog of the Duhamel formula (fl8|) and (|2T|) for the second-order cumulant 
2t2 ( — t, Yi, s + 1), we have 

t 



i r 5+1 

lim||(-2l 2 (-t,n,s + l) - / d* 1 IJ^i(-* + *i>i) x 

o i =1 



s+l 

i=l Z=l 



( - ^au*, * + 1)) n&(-*i, o)/-+iL (w . +1) = o. 



s 

lim||(^2l 1+n (-t) - /"dti... / ^nnSiH + *i,j)^(-A/Lt(*i,s + l)) 

^ 6 n n .7=1 il=l 



In general case the following equality holds: 

t tn—l 

& S 

" 1 ' ^ (-AT int (2i,S + 1)) X 

b o ■ J=1 il= " 

s+l s+n— 1 s+n— 1 

X J| +t 2 ,jl) ••• ] [ 0l(-*n-l+*»,Jn-l) Yl (~ Mnt(in,S + n)) X 

il=l J»-l=l *n=l 

x II^i(-*» > i»))/* + «L CK ^ ) = 0. (28) 

Thus, if, for the initial data F s (0) G £ 1 ('H S ), there exists the limit / s (0) G £ 1 (7Y S ), i.e., 

lim|| e ^ s (0)-/ s (0)|| £l(Ws) = 0, 

then, according to (|28| for an arbitrary finite time interval, there exists the mean-field limit 
of solution ([TBI of the BBGKY hierarchy 

hm\\e s F s (t)-f s (t)\\ mHs) = 0, 

where f s {t) is given by the series 

t t n -i 

oo „ „ s 

f s (t,l,...,s) = ^ / dti... / d£„Tr s+1) ... )S+n JJ^i(-t + ti, j) x 

n=0 ^ ^ ,=1 

s s+l s+n— 1 

x ( --Mnt(«i,s + 1)) JJ Gi(-h +t 2 ,ji) ■ ■ ■ Yl Gl( 

H=l ii=l in-i=l 

s+n— 1 s+n 

x (-^(4, s + n)) Y[ £i(-t„,j„)/s+n(0), (29) 
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which converges for a bounded interaction potential on a finite time interval [111 ]. 

If /(0) G £q(JF w ), the sequence /(i) = (/, fi(t), . . . , / s (t), . . .) of limit marginal density 
operators (I2"$j) is a strong solution of the Cauchy problem of the Vlasov hierarchy 

p) s s 

g~ t fs(t) = E ( " Na(i))f s (t) + J] Tr s+1 ( - AT int (z, S + l))/ s+1 (t), (30) 

i=l i=l 

fs(t)\t=0 = W), 8>1. (31) 

We observe that, if the initial data satisfy the chaos property (for particles obeying 
Maxwell-Boltzmann statistics) 

s 

f a (t,i,...,s)\ t=0 = l[f 1 {o,j), s>2, 

then solution ( 1291) of the initial- value problem of the Vlasov hierarchy ( I30l) - (13T1) possesses of 
the same property 

s 

f.(t,l,...,8) = l[f 1 (t,j), s>2. (32) 

3=1 



To established equality (I3"!2"j) . we introduce marginal correlation density operators [311 ] 



G.(t, 1, . . . , s) = ^ Tr s+i,...,s + n a a+n (-t, 1, . . . , s + n) G x (0, i), (33) 

n=0 ' i=l 

where 2t s+n (— £) = 2l s+ri (— £, 1, . . . , s + n) is the (s + n)th-order cumulant flT5l) of the groups 
of operators (Tl3l . and Gi(0) = i*i(0). In the same way as (|28l) for arbitrary i 6 R, we 
establish the equality 

S3ll^ 2t «+n(- t ' 1 '---' S + 7l )/-+»L(W. +B ) = - ( 34 ) 

Let 

hold. Then, according to (IMj) for the correlation density operators ( 1331) . we obtain 

!™lk G .(*)L<*) = - ( 35 ) 
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In view of the fact that the marginal density operators ffl6|) are expressed in terms of the 
correlation density operators fl33|) by the cluster expansions 

s 

F s (t 7 Y) = ]jF l (t,i)+ HG\Xi\&Xi)> s > 2 > 

1=1 P : {Y} = \J i X i , X ' CP 

|P| ^ s 

and taking equality (I35p into account, the following statement is valid. 
If there exists the mean-field limit of the initial data F s (0) G £}(H S ) 

s 

Km \\e s F s (0, l,..., a )-n = 0, 

3=1 

then, for a finite time interval for solution l[Tb)) of the BBGKY hierarchy the limit 

s 

Km \\e s Fs(t, 1, , s) - H fi(t,j)\\ sl{Ht) = 

3=1 

holds, where fi(t) is the solution of the Cauchy problem of the quantum Vlasov equation 

j^/i(M)= (-Ar (l))/ 1 (t,l) + Tr 2 (-A^ nt (l,2))/ 1 (t,l)/ 1 (t,2), (36) 

/i(t)|t=o = /i(0). (37) 

Thus, in consequence of the chaos property ( 1321) . we derive the quantum Vlasov kinetic 
equation (1361) . 



B. On the nonlinear Schrodinger equation 

For a system in the pure state, i.e. fi(t) = \ipt)(^>t\ ( Pi/> t = \i ; t)('4 , t\ is a one-dimensional 
projector onto a unit vector \ipt)) or in terms of the kernel f\(t, q, q') = ip(t, q)i/j(t, q') of the 
marginal one-particle density operator fi(t), the Vlasov kinetic equation (1361) is transformed 
to the Hartree equation 

g) = --A^(t, q) + / <V$(g - g'Mit, q') | 2 ^(t, g ). (38) 



dt Ty ^' 2 

If the kernel of the interaction potential $(g) = 5(g) is the Dirac measure, then from ([38 
we derive the cubic nonlinear Schrodinger equation 

i^(t,q) = -^A y ^(t,q) + \ifj(t,q)\ 2 ^(t,q). 
18 



Thus, the following statement holds: 

lim||e s F s (t)-|^)(^r s || £1(Hs) = 0, 

where \ipt) is the solution of the cubic nonlinear Schrodinger equation. 

In the case of representation ffl~9]) of the solution of the Cauchy problem (fTT])-(|12D of the 
BBGKY hierarchy by the iteration series the last statement is proved in works 

C. The mean- field limit of the generalized quantum kinetic equation 



We construct the mean-field limit of a solution of the initial-value problem of the gener- 
alized kinetic equation (|25|) . 

If there exists the limit /i(0) e £}{T-L\) of initial data 



lim || C F X (0) -/i(0)|| £l(Wi) = 0, 

then, according to ( 1271) and ( 1281) for an arbitrary finite time interval there exists the limit of 
solution (1241) of the generalized kinetic equation (1251) 



\mv\\eF x {t)- f x {t)\\ zl(:Hi) = ^ (39) 



where f\{t) is a strong solution of the Cauchy problem (I3"oT) - (13"TI) of the quantum Vlasov 
equation represented in the form of the expansion 

t *n— 1 

00 ™ ™ s 

f 1 (t,l) = ^2 / dtx... / dt n Tr s+1 _ s+n Y[g 1 (-t + t 1 ,j)x 

n=0 { { j=l 

s s+1 s+n— 1 

x^(-M Bt (i 1 ,8 + l))Y[Q 1 (-t 1 + t 2 J 1 )... Yl Gi(-t n -i + t n ,j n _i)x 

ii=l ji=l j n -i=l 

s+n s+n 

x ( - Mtoiin, s + n)) JJ Qx{-t n ,j n ) n h(0, i), (40) 

i„=l jn = l 1=1 

and the operator A/lnt is defined by formula (j3j). For bounded interaction potentials, series 
(|40p converges for a finite time interval. 

If f s G i2 1 (7^ s ) and an interaction potential is a bounded operator, then for scattering 
operators (123]) an analog of the Duhamel formula holds: 

Is s 

(G s (t,l,...,s)-I)f s = e / drl[g l (T,l)(- Mnt(i,j))Gs(-r)f s . (41) 

o 1=1 *<i =1 
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Then, according to definition (12T1) of the evolution operators 2Ji +n (t, {1, . . . , s}i, s+1, . . . ,s+ 
n), n > 0, from expansion (120]) and equality (j4"Tl) we establish 

Iim||(fn 1 (t,{l,...,a} 1 )-/)/ - |L fW .=0. 



Correspondingly for n > 1, it holds: 



lim||9J 1+n (t)/ s+n || £l(Ws+ji) = 0. 



Since a solution of the initial- value problem (I231)- (l26l) of the generalized kinetic equation 
converges to a solution of the initial- value problem (!36|) - (l37|) of the quantum Vlasov kinetic 
equation as (1391) . for functionals (1201) we obtain 

s 

Jim ||e s F s (t, 1, . . . , s | F 1 (t)) - J] A(t, j)|| £l(Ws) = 0, 

i=i 

where is defined by series ( 140|) which converges for a finite time interval. The last 

equality means that the chaos property fl32l) preserves in time in the mean-field scaling 
limit. 

Thus, we conclude that the results of the previous subsection concerning the derivation 
of the Hartree equation and the nonlinear Schrodinger equation take place also in the case 
of the generalized quantum kinetic equation (1231) . 

D. The mean-field limit of the dual BBGKY hierarchy solution 

Consider the mean-field limit of a solution of the initial-value problem of dual BBGKY 
hierarchy (J2J). 

For an arbitrary finite time interval, there exists the following limit of the *-weak contin- 
uous group of operators (j7j) in the sense of the *-weak convergence of the space £(Ti s ) 

s 

w*-lim (g s (t)g s - Y[Qi{-t,j)g s ) = 0. (42) 

3=1 

According to an analog of the Duhamel formula (1T8|) and (|42|) for the second-order cumulant 
2l 2 (t, 1, 2) in the same sense as above, it holds: 

* 2 2 

w*-lim(-2t 2 (t,l,% 2 - I dt 1 f[g 1 {t-t 1 ,j)Mnt{l,2)f[g 1 {t 1 ,l)g 2 ) =0. (43) 

{ 3=1 1=1 
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Thus, if, for initial data G s (0) G £(H S ), there exists the limit g s (0) G £,(H S ), i.e. if it 
holds: 

w*- lim (e- s G s (0) - g s (0)) = 0, (44) 

then, according to fj42l and (j4"3]) for an arbitrary finite time interval there exists the mean- 
field limit of solution (TlOl) of the dual BBGKY hierarchy ([2]) in the sense of the *-weak 
convergence of the space £(TC S ): 

w*-hm(e- s G s (t)-g s (t))=0. (45) 

The limit operator g s {t) in (j4"5]) is given by the expansion 

g s (t,Y) = J2 I dU... j dt n g° s (t - h) Mnt^^k^GtAti-h)... 

s 

)gl n (t n )g s ^ n (0,Y\{i kl ,...,i kn }), (46) 

where 

Gs-n+lttn-l - = £°- n +l(*n-l - 4, ^Yfe, • • • , k n -i}) = Gl(t n -1 ~ t n , j) 

jGy\{j fel ,...,j fcn _ 1 } 

is the group of operators ([7]) of noninteracting particles. 

If g(0) G £(JF W ), the sequence g{t) = (g , g±(t), . . . , g s (t), . . .) of limit marginal observables 
( 1461) is a generalized solution of the initial-value problem of the dual Vlasov hierarchy 

ft s s 

-g s (t,Y) = J2^9s(t,Y)+ Mnt(ju32)g s -i(t,Y\{j 1 }), (47) 

g s (t)\t=o=9s(0), s>l. (48) 
Consider the mean-field limit of the additive-type observables, i.e. 

G«(0) = (0,gS 1} (0,1),0,...). 
In that case, solution fflUj) of the dual BBGKY hierarchy (T5]) has the form 

s 

G^(t,Y) = ^ s (t,Y)Y,G?(0,j). (49) 

3=1 
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If, for the additive-type observables G^(0), condition (Jill) holds, i.e. 

w*-lim {e-^W-gPip)) =0, 
then, according to statement (1451) . for ( 1491) we have 

w*-hm( e - s G«(t)-^(t))=0, 

where 

g?(t,l) = g 1 {t,l)g?\0,l) i 

I 2 2 







3=1 1=1 



or as a special case of ( 1461) the limit operator gi^(t) is defined by the expansion 

t t s -2 s 

gi 1 \t,Y) = J dt,... J dt s ^g° s {t-t x ) Yl Mnt(%,4 2 )d(tl-t2)... 

*fc 1 7^*fc 2 = 1 

• • • £ 2 °(t s -2 - Mntiik^, ih,)Qi{t-i) 9i\0, Y\{iki, ■ ■ ■ , (50) 



Let the initial state satisfy the chaos property ( 1321) 

s 

/i e) (0,l,...,a) = n/i(0,j), s>2. 

i=l 

Then, if g(t) G £ 7 (JF W ) and /i(0) G £ 1 (Hi), the mean value functional 

oo 1 s 

(g(t), f(0)) = ~ B 1, , a) ]J /i(0, i) 



s=0 i=l 



exists, provided that U/i^H^u?^) < 7. 
In consequence of the equality 



00 1 s 

{g {1) (t), j< c \o)) = £ - Tri,..., s gS-Ht) J] /i(o, = Tn ^(o)/^, 1), 



S! 

s=0 i=l 



where <7^(£) is given by ( 1501) and 1) is solution ( )40|) of the quantum Vlasov equation 
( 1361) . we find that the initial- value problem (l4"T]) - (l4"8]) for additive-type observables and the 
initial state /^(0) describes the evolution of quantum many-particle systems as by the 
Vlasov kinetic equation. 



22 



Correspondingly, the chaos property (1321) in the Heisenberg picture of evolution of quan- 
tum many-particle systems is fulfilled, which follows from the equality (k > 2) 

{g {k) (t), / (c) (o)) = £ - n A(o. = ^ Tr i,..,^f(o) n A(*. 0, 

s=0 ' i=l ' ' i=l 

where fi(t, 1) is given by expansion (HOl) . 

Thus, if the initial state is a pure state, i.e. / s (0) = \i(Jo) {ijJo\® s , we conclude that in 
the Heisenberg picture of evolution, the initial-value problem ( )47l) - (l48l) of the dual Vlasov 
hierarchy describes the evolution of quantum many-particle systems which is governed by 
the Hartree equation ( 1381) in the Schrodinger picture of evolution or it is governed by the 
cubic nonlinear Schrodinger equation, if the interaction potential = 5(q) is the Dirac 
measure. 



IV. CONCLUSION 

The concept of cumulants (Q of the groups of operators (j7]) of the Heisenberg equations 
or cumulants (fl5l) of groups of operators (fl3l) of the von Neumann equations forms the basis 
of the groups of operators for quantum evolution equations as well as the quantum dual 



BBGKY hierarchy and the BBGKY hierarchy for marginal density operators 28(], 29j, 32]. 

As was mention above, for the initial data F(0) e £^(JF W ), the average number fTTTj) of 
particles is finite. In order to describe the evolution of infinitely many particles [5] , we have to 
construct solutions for initial marginal density operators belonging to more general Banach 
spaces than £^,(JF W ). For example, it can be the space of sequences of bounded operators 
containing the equilibrium states |4|. In that case, every term of solution expansions for 
BBGKY hierarchy ( iTTjl and correspondingly for the generalized quantum kinetic equation 



,[8 



In the case of the dual BBGKY 



( 125]) and functionals (l20j) contains the divergent traces [5j 
hierarchy (j2J), the problem consists in the definition of mean value functional (J6]) which 

26(,[32|. The analysis of such a question for quantum systems 



contains the divergent traces 
remains an open problem. 

We have formulated two new approaches to the rigorous derivation of quantum kinetic 
equations from underlying many-particle dynamics. These approaches enable one to describe 
the kinetic evolution if the chaos property (|32j) is not fulfilled initially, i.e. in the presence 
of initial correlations. Such Cauchy problem takes place in the case of kinetic evolution of 
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the Bose condensate . As a result of these approaches, we can formulate the kinetic 
equations both for a Bose gas and a Bose condensate, i.e. the nonlinear Schrddinger equation 
and the Gross-Pit aevskii equation [ill]. [121]. respectively. 
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